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ABSTRACT 

In this paper, we address how block rotation and complex distributed de- , 
formation in the Earth’s shallow crust may be explained within a stationary 
regional stress field. Distributed deformation is characterized by domains of 
sub-parallel fault-bounded blocks. In response to the contemporaneous activity 
of neighboring domains some domains rotate, as suggested by both structural 
and paleomagnetic evidence. 

Rotations within domains are achieved through the contemporaneous slip 
and rotation of the faults and of the blocks they bound. Thus, in regions of dis- 
tributed deformation, faults must remain active in spite of their poor orientation 
in the stress field. Traditional friction models cannot account for this mecha- 
nism. To solve this problem we developed a model that tracks the orientation of 
blocks and their bounding faults during rotation in a 3D stress field. Mechani- 
cally, we considered Coulomb criteria for rock fracture, as an upper bound, and 
fault slippage, as a lower bound, between which block rotation is expected. 

In our model, the effective stress magnitudes of the principal stresses (< 7 j, <Ji , 

03 ) are controlled by the orientation of fault sets in each domain. Therefore, ( 1 ) 
adjacent fault sets with differing orientations may be active and may display dif- 
fering faulting styles, and ( 2 ) a given set of faults may change its style of motion 
as it rotates within a stationary stress regime. The style of faulting predicted 
by our model depends on a dimensionless parameter <f> = (< t 2 — cr 3 ) /(< 7 j — <73). 
Thus, we present a model for complex distributed deformation and complex 
offset history requiring neither geographical nor temporal changes in the stress 
regime. 

We apply the model to the Western Transverse Range domain of Southern 
California. There, it is mechanically feasible for blocks and faults to have ex- 
perienced up to 75° of clockwise rotation in a <f> = 0.1 strike-slip stress regime. 
The results of our model suggest that this domain may first have accommo- 
dated deformation along preexisting NNE-SSW faults, reactivated as normal 
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faults. After rotation, these same faults became strike-slip in nature. Subse- 
quent rotations could have resulted in the present day E-W high angle reverse 
faults. This history agrees with both prominent structural phases documented 
for post-Oligocene activity and paleomagnetically inferred rotations of this do- 
main. 
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CRUSTAL DEFORMATION BY BLOCK ROTATION 

1. The problem 

Distributed crustal deformation is characterized by sub-parallel sets of faults 
(Freund, 1970; Freund, 1971; Garfunkel, 1974; Luyendyk et al., 1980; Ron et al., 
1984; and others). Sets are distinguished from their neighbors by orientation and 
often by faulting style as well. In a region like Southern California, throughgoing 
faults often mark the boundaries of these domains such as the San Andreas and 
Garlock fault shown in Fig. 1. 

If the regional stress field is homogeneous and stationary throughout a region 
like Southern California (Zoback et al., 1987), then how can we have active fault 
sets in such varying orientations? At least three solutions can be suggested: 

(1) In regions of distributed deformation, fault sets behave in accordance with 
unknown friction criteria. 

(2) The stress field is not homogeneous and stationary. Instead, it changes 
orientation from one domain to the next. 

(3) The stress field is homogeneous and stationary and the fault sets slip in 
accordance with friction criteria but in some domains blocks and faults rotate. 

The answer most probably lies in a combination of these three extreme cases. 
Fault behaviour is most likely a function of the slip rate and the history of 
slip; the stress field must be inhomogeneous to some extent across a region of 
distributed deformation; rotation of blocks and faults must indeed occur when 
many domains of fault sets coexist. 

In this paper, we restrict our attention to the third case: the rotation of 
faults and of the blocks between them. Accumulating paleomagnetic and struc- 
tural evidence shows that in region of distributed deformation, many structural 
domains have rotated in the past, and some are rotating today (see Nur et al., 
1986; Ron et al., 1988). Different block rotation mechanisms have been proposed 
(see Molnar, 1988). Here, we propose a 3D block rotation mechanism based on 
Nur et al. (1986) 2D model. For this reason, we first give a brief summary of 
the kinematics and mechanics of the 2D version. 

2. Kinematics of block rotation in 2D 

The kinematics of block rotation in strike-slip tectonic regimes, were orig- 
inally proposed by Freund (1970; 1974) on the basis of structural data alone. 
He observed two sets of strike slip faults with an angle of 145° between them 
measured in the maximun compressional direction. He attributed this angular 
spread to the tendency of strike slip faults to rotate about a vertical axis away 
from the maximum compressional direction. In fact, rotation of faults as a mech- 
anism for accommodating deformation, was first recognized in normal tectonic 
regimes (Ransome et al., 1910; see references in Jackson and White, 1989). 

This rotation is illustrated schematically in Fig. 2. In this simple case, the 
vertical stress S v = = gravity , is the only stress acting. As the books slide on 

the shelf they rotate away from Uj, the direction of maximum compression: in 
a counter-clockwise (CCW) sense (domain A) or clockwise (CW) sense (domain 
B) depending on the orientation of the books. From this example, one can see 
that — whether strike-slip or dip-slip — right-lateral fault motion leads to CCW 
sense of rotation, and left-lateral movement leads to CW sense of rotation of 
blocks and faults. 

The block rotation model assumes that the rotating blocks are rigid. Luyendyk 
et al. (1985), Ron et al. (1984), Carter et al. (1987) and Terres and Luyendyk 
(1985) have validated this assumption. In the domains they studied, two sig- 
nificant correlations were found. The first one between the measured sense of 
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slip, the expected sense of rotation and the paleomagnetically inferred axis of 
rotation (as explained in Fig, 2); the second one between the known amount 
of total displacement across a fault set, the measured average spacing between 
sub-parallel faults, and the amount of rotation inferred paleomagnetically for a 
domain. 

3. Mechanics of block rotation in 2D 

In Freund’s model, as faults slip, the blocks they bound rotate away from 
the direction of maximum compression. Nur et al. (19861 added a mechanical 
limit to this kinematic model for the case of vertical striKe slip faults rotating 
in a strike-slip regime. In Nur’s block rotation (BR) model, two additional 
constraints are introduced. One is given by the Coulomb criterion for sliding, 
and the second one by the Coulomb criterion for fracturing (eq. (1) and (2) in 
the Appendix). 

Within a region of distributed deformation, where a single fault set orienta- 
tion characterizes each domain, the first constraint implies that fault sets may 
remain active even if poorly oriented in the stress field, and the second con- 
straint sets a limit on how poorly oriented the fault set may become. A fault 
is considered optimally oriented when the intermediate stress a 2 is contained in 
the plane of the fault and the shear stress required for slip along the fault is 
such that the crj — cr 3 Mohr circle is just tangent to the sliding fine (Fig. 3). 

In Nur’s model, faults rotate within a stationary stress field. This implies 
that during rotation the stress magnitudes must change in a domain, for blocks 
and faults to continue to slip and rotate (Fig. 4). When the magnitude of the 
differential stress reaches the strength of the intact crust (the fracture line on the 
a n — r plot), a new set of more optimally oriented faults forms, and the old set 
becomes locked. As a result, a discrete range of fault orientations is predicted 
by the BR model. Only when sufficiently large rotations occur in a domain will 
cross-cutting generations of fault sets give rise to a complex pattern of faulting, 
such as those observed in situ (Angelier et al., 1985; Ron et al, in press). 

The mechanical constraints of the BR model provide an important step to- 
ward quantitative estimates of the contribution of block rotation to distributed 
deformation in the Earth’s crust. The relationship among friction, strength and 
the amount of rotation that a single set of faults can experience are detailed in 
Nur et al. (1986). 

As Fig. 4 shows, block rotation can cause poor fault orientation. However, 
the 2D model is limited in its application since the intermediate principal stress 
<r 2 must remain in the fault plane. This limitation prevents any change of faulting 
style during rotation. A 3D formulation is necessary to model the more common 
case of rotation along reactivated faults, where a 2 is not necessarily in the plane 
of the fault. Then, we can identify the conditions that allow a fault set to 
drastically change its behaviour as it rotates. 

4. Block rotation in 3D 

The focus of this paper is to understand distributed deformation and the 
relationship between fault slip and the rotation of faults and blocks in a three 
dimensional stress field. The fundamental strenght and friction criteria proposed 
in the 2D model are preserved. In addition, the 3D model considers all faulting 
styles that can occur in a general stress field for any starting orientation of 
modeled faults. 

We will first discuss the more fundamental assumptions that must be made 
when modeling distributed deformation and block rotation in a 3D stress field. 
Then, we will discuss the implications of combining the 3D stress Mohr circle 
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with friction criteria applied to the present-day tectonics of Southern California. 
In the main part of the paper, we will discuss the principles of the 3 D BR model. 
As the results will show, even a simple 3D BR model can produce complicated 
faulting histories and complex distributed deformation geometries. 

Because faults generally change their behaviour as they rotate in a station- 
ary stress field, a model of faulting and rotation in 3D is fundamental to our 
understanding of this mechanism. In the final part of the paper, as a practi- 
cal application of our model, we will compare the known rotation history of the 
Western Transverse Range domain in Southern California with the results of our 
model. We find encouraging similarities between our 3 D BR model predictions 
and the complex structural and geological record of this domain. 

THE STATE OF STRESS 

1 . One principal stress assumed vertical 

Following Anderson (1951), the vertical stress S v is assumed to be a principal 
stress. The magnitude of S v determines the tectonic stress regime: normal if 
S v — <7i, reverse if S v — o 3 , and strike-slip if S v = 02 — where 0 \ > <Ji > a 3 and 
compressive stresses are positive. Note that beyond Anderson’s view, a given 
tectonic stress regime does not imply a specific faulting style. Particularly in 
regions of distributed deformation because they are characterized by domains of 
fault sets. 

2. Stress models 

As shown in Fig. 4 for the 2D case, stress magnitudes must change in a 
domain because we assume that faults and blocks rotate within a stationary 
stress field. Since we do not know how this change may take place, we must 
make a number of simplifying assumptions. One possible assumption is that the 
volume of the crust being deformed remains constant so that hydrostatic stress 
is constant 


<7i + cr 2 + cr 3 = constant 


Another possibility is that the stress ratio (f> remains unchanged 


<f> = 


<72 ~~ <7 3 

— constant 

0\ — O3 


The stress ratio <f> = (&2~ cr 3)/{ (7 i —03) is often used as a dimensionless parameter 
to describe the 3D state of stress in the Earth’s crust. It may vary from <j> = 0 
when a 2 = o 3 to <j> — 1 when a 2 = 0\. 

3. Fault geometry representation in 3D 

Whatever assumptions are made, we need a way of tracking fault plane orien- 
tations during rotation in a 3D space. In this paper, we use two representations: 
the Mohr circle, which plots the orientation of faults in a o n — t plane, and 
the Wulff-projection stereonet, that plots the orientation of fault planes in the 
principal stress axis reference. 

THE 3D MOHR CIRCLE 

In 2 D space, a single angle defines fault orientation. In 3 D space, two angles 
are required to define a pole to a fault plane. Fig. 5 a shows how these angles are 
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represented in the Mohr circle. For a more complete treatment of the 3D Mohr 
circle representation we refer the reader to Jaeger and Cook (Chp 2.6, 1969). 
Let us consider here only three extreme cases: 

(1) pole P 2 which falls exactly on the U\ — <73 circle: this is equivalent to the 2D 
representation in which <72 is in the plane of the fault. 

(2) pole Pi which falls exactly on the a 3 - a 2 circle: in this case cr x is in the 
plane of the fault and 

(3) pole P 3 which falls exactly on the a 2 — o x circle: in this case er 3 is in the 
plane of the fault. 

When a pole (PI falls in the shaded region in Fig. 5a, then all three principal 
stresses are off the plane of the fault. In this case ori and a 3 , the angles that the 
fault plane normal makes with a x and cr 3 respectively, are calculated by drawing 
two circles concentric with the two small o x — <72 and <7 2 — <7 3 circles and passing 
through the pole. Then the intersection of the 0 \ — <t 2 concentric expansion 
with the a 2 — a 3 circle defines a 3 , and the intersection of the a 2 — <7 3 concentric 
expansion with the — a 2 circle defines a x as shown in Fig. 5a. The Mohr 
circle representation is powerful because it represents stress magnitudes, friction 
criteria and the geometry of faults in the stress field, all in one graph. 

THE WULFF-PROJECTION 

The same fault normals shown in Fig. 5a axe plotted on a Wulff lower 
hemisphere projection in Fig. 5b: P 2 plots on the line joining the a 1 — cr 3 
directions, P x plots on the line joining the a 2 — <j 3 directions and P 3 plots on the 
line joining the cr 1 — a 2 directions. Depending on the stress regime, either <7 1? o 2 
or a 3 is represented in the down orientation throughout the paper. This type of 
representation is often used by structural geologists because it makes it easier 
to visualize the geometry of faults in the stress field. 

4. Choice of material parameters 

Finally, we need to assume some values for the mechanical constraints of the 
BR model. The relationship between material parameters and block rotation is 
discussed in detail in Nur et al. (1986). As we will see for the case of the Western 
Transverse Range, the specific values do not affect the qualitative aspect of the 
results. In this paper, we are more concerned with the changing behaviour of 
faults as they rotate. In the model, we assume the following values for the 
material parameters (Handin. 1969): 1.0 for the coefficient of friction of intact 
rock (it could be much greater) and 0.6 for the coefficient of friction of preexisting 
faults fit could be less). 1000 bars for the cohesion of intact rock (could be 
greater) and 50 bars for the cohesion of preexisting faults (could be 0). Given 
these values, the stress limits (Coulomb criteria) employed in the model become: 

(1) The upper limit, representing the strength of the crust (the fracture line) 

r 0 = 1000 + 1.0<7„ (in bars ) 


(2) The lower limit, representing the strength of preexisting faults (the sliding 
line) 

tj = 50 + 0.6<7 n (in bars ) 

(3) the tensile limit, assuming the crust cannot withstand any tensile stress 

cr 3 > 0 (<7 3 < a 2 < <7i) 

3D MOHR CIRCLE: A KEY TO MIXED STYLES OF FAULTING 
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In a three dimensional situation, the value of the intermediate stress plays 
a key role in determining the style of faulting. This contrast with the two di- 
mensional one where <7 2 is always in the fault plane. By itself, this result is 
not new ( Bott, 1959; McKenzie, 1969; and others). However, here we include 
cr 2 specifically to allow us to understand how active faults become poorly ori- 
ented during 3D block rotation. Before considering block and fault rotation, 
let us first enlist the help of the Mohr circle to see when preexisting faults can 
be reactivated — and what this means for regionally distributed deformation in 
Southern California. 

1. Friction criteria and the 3D Mohr circle 

Clearly, reactivation of preexisting faults is constrained by the strength of the 
crust, the strength of the faults and the relative magnitudes of the three principal 
stresses. Jaeger and Rosengren (1969) discuss the influence of these parameters 
in more detail. Let us consider the material parameters values discussed in the 
previous section and a a 3 /c r x ratio exceeding the sliding line but not the fracture 
line, as shown in Fig 6. 

Consider only three sets of preexisting faults as sketched in Fig. 6: (a) set 1 
which contains the <7i direction, (b) set 2 which contains the <7 2 direction and (c) 
set 3 which contains the <r 3 direction. The faulting styles expected along each 
fault set are also sketched. They depend on the tectonic stress regime and the 
orientation of each fault set. The normals to actively slipping faults plot within 
the shaded region of the Mohr circle. The size of the shaded region depends on 
the 03/(7 1 ratio and the <f> value. Assuming low o 3 j<J\ ratio, at lower 4> values, 
say less than 0.4, all three fault sets will be active. In the normal regime, set 1 
will be reactivated as strike-slip faults, set 2 and set 3 as normal faults. In the 
reverse regime set 1 and set 2 will be reactivated as reverse faults and set 3 as 
strike-slip faults. In the strike-slip regime, set 1 will be reactivated as normal 
faults, set 2 as strike-slip faults and set 3 as reverse faults. At <f values exceeding 
0.4, fault set 3, which contains <73 in its plane, will be locked, while the other 
two fault sets can be reactivated. 

These cases represent the most extreme styles of faulting that can be expected 
in the three tectonic stress regimes. Clearly in a domain characterized by an 
obliquely slipping fault set, movement could occur if it plotted within the shaded 
region of the 3D Mohr circle. A summary table for the limiting cases is provided 
in Table 1. 


Stress Ratio 

Tectonic stress regime 

Faults 

normal 
S v = (T\ 

reverse 

S v = <73 

strike-slip 

Sy = & 2 

low 

normal 

reverse 

strike-slip 

mmm 

0 . < <j> < 0.4 

strike-slip 

reverse 

normal 

mm 


normal 

strike-slip 

reverse 

set 3 

high 

normal 

reverse 

strike-slip 


0.4 < <f> < 1.0 

strike-slip 

reverse 

normal 

mm 


Table 1. Limiting cases of faulting styles expected in three tectonic stress regimes 
as functions of the stress ratio <j) = (<7 2 — o 3 )j(cr\ — 03 ) — for each of the three 
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sets aligned with one of the principal stresses in their planes and for a low as/a-i 
ratio. 

2. Southern California domains: an example of distributed deformation 

Southern California is characterized by domains of faults separated by through- 
going boundary faults (Fig. 1). The shaded regions in Fig. 1 represent domains 
where block rotation has been inferred (Luyendyk et al., 1985) and therefore 
where we expect faults to be poorly oriented in the present-day stress field. 
Whatever the stress field orientation, we do observe a mixture of faulting styles. 
In the Western Transverse Range domain (WTR), high-angle reverse-oblique 
slip along E-W trending faults has been documented (Lee et al., 1979; Yerkes 
and Lee, 1979a, b). In the Mojave domain (MOJ), Sauber et al. (1986) have 
described right-lateral strike-slip along NW-SE trending vertical faults. In the 
East Transverse Range domain (ETR), activity is left-lateral strike-slip along 
E-W trending vertical faults (Jones, 1988; Powell, 1982). 

Can this diverse fault behaviour be the result of a regionally stationary stress 
field — in accordance with friction criteria? Note first how the San Andreas strike- 
slip system dominates the region’s tectonics. Consider next, from the arguments 
above and Table 1, that reverse and strike-slip faulting may coexist within a low 
<f) strike-slip stress regime and low a z /a x ratio. Lastly, we need an estimate of 
the < 7 1 direction and of the value of /i/, the coefficient of frictional sliding, to 
apply the 3D BR model to this region. Both, the cr 1 direction and the fij value, 
are constrained by the orientation of strike-slip faults in the MOJ domain and 
of high angle reverse faults in the WTR domain: a x must trend N20°E and 
(*/ = 0.4 for MOJ and WTR faults to slip. Representative fault plane solutions 
shown in Fig. 1 are plotted on the 3D Mohr circle of Fig. 7. The plot is 
normalized by the value of cr x , assuming a depth of faulting of 5-10 kms, an 
orientation of a x trending iV20 6 E, a \i / = 0.4 and a 4> = 0.1. Faults of the ETR, 
as well as those of the WTR and MOJ domains, plot closely to the sliding line. 
Therefore, in our simple model of distributed deformation, faults sets in all three 
domains can slide in accordance with friction criteria in a regionally stationary 
strike-slip regime. 

The BR model can not presently account for slip along domain bo un ding 
faults (Fig. 1). In fact, the model predicts that very little shear stress is resolved 
on the San Andreas fault, the major throughgoing boundary fault in this region. 
While this agrees with recent borehole findings at Cajon Pass (Zoback et al., 
1987), questions remain concerning the mechanics of lithospheric faults, but they 
are beyond the scope of this study. 

In summary, we see discrete fault orientations that separate into domains. 
In Southern California, even poorly oriented fault sets can slip in a a x jV20 °E 
directed strike-slip regime with <f> = 0.1 and fij = 0.4. We will return to these 
values in the final part of the paper when discussing the faulting history of the 
WTR domain. 

A fundamental question must be considered at this point: how do faults 
that presumably form in a favorable orientation become poorly oriented? As 
mentioned previously, extensive paleomagnetic studies (references in Luyendyk 
et al., 1985 and Dokka, 1989) indicate that these domains have undergone a 
complex history of deformation characterized by rotations about a vertical axis. 

In the following section, we hope to demonstrate that block rotation in a 3D 
stationary stress field provides a simple mechanism to accommodate distributed 
deformation that explains, in accordance with friction criteria, those rotations 
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that have been documented. 

3D BLOCK ROTATION 

In the previous section, we defined conditions in a 3D stress field under 
which reactivation of poorly oriented faults can occur (Table 1). We showed 
that domains of poorly oriented faults can be active in a regionally stationary 
stress field if the limits for the differential stresses in the Earth’s crust lie between 
the upper bound (rock strength) and the lower bound (fault strength). 

We will now investigate a process which allows faults, to pass from optimal 
to poor orientation; this is the rotation of blocks and their associated faults. 
Block rotation was modeled kinematically in 2D by Freund (1974), Garfunkel 
(1974) and Ron et al. (1984) and mechanically in 2D by Nur et al (1986). In 
this paper we present a 3D mechanical modeling program. Given strength and 
friction values, our model computes both the maximum rotation possible for a 
given fault set, and the conditions needed to induce a change of faulting style 
for that set. 

1. The model 

The block rotation model (BR) describes idealized domains of fault sets in 
situ. Faults in a set are typically sub-parallel, so it is reasonable to consider 
the rotation history of one fault to be representative of the entire set. In the 
BR model, we assume that all faults in a given set are active simultaneously. 
Thus, deformation remains uniform throughout the domain. We also assume 
that the rotating blocks are rigid. Faulting is therefore brittle. We expect a a 
detachment at depth (ca. 15 kms), that decouples upper crustal rotations from 
lower crustal ductile shear. 

In the brittle upper crust, we require four main assumptions for the 3D BR 
model. The first, a mechanical constraint, assumes that the Coulomb criterion 
controls sliding of faults in the set. The second assumes that fault slip is directed 
along the maximum resolved shear stress for the fault. The third, a kinematic 
constraint, assumes that both blocks and fault planes rotate away from the <j\ 
direction. The fourth and final assumption is the stationarity of the principal 
stress directions. Thus, once fault sets have slipped and blocks rotated, stress 
magnitudes must change to allow further slip and rotation. 

To facilitate the analysis, we assume that the magnitude of the principal 
vertical stress S v remains constant throughout rotation. The other two stress 
magnitudes vary to conserve stress ratio <f> (see Appendix for discussion). As 
mentioned earlier, more general stress histories are possible but they require 
additional information. Here, we present results for a constant <f> stress history 
only. 

Given a stress model and a stationary direction of the stress field, we must 
keep track of the fault orientation during rotation. Two new concepts are in- 
troduced: the stress path and the rotation path of a rotating fault set. These 
concepts form the basis of the 3D BR model. Estevez et al. (1987) discussed 
the stress path concept in considerable detail. The derivation is summarized in 
the Appendix here. This paper will focus mainly on the rotation path concept. 
The rotation path of a rotating fault set can be represented on a Mohr circle, 
as the successive (<r n , r) values of the fault set on the Coulomb sliding line, and 
on the Wulf-projection, as the subsequent positions of the fault set in the stress 
field. 

2. 2D rotation path 

The rotation path, a fundamental concept of 3D fault rotation, is illustrated 
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schematically on a 2 D Mohr circle in Fig. 8 . Assume that a given domain 
contains one preexisting fault set with its normal aligned at c*i from the <J\ 
direction (Fig. 8 a, point 1 ). For the fault set to slip, the stress must rise from 
point 1 to point 2 . By the Coulomb criterion, shear stress at point 2 can drive 
the fault set into motion. At slip, the fault set rotates, eventually assuming 
an orientation «2 (Fig. 8 b). The angles c*i and a 2 represent fault normals! 
Therefore, as blocks and faults rotate away from < 7 j, normals to faults approach 
the <7 1 point along the a n axis — as shown earlier in Fig. 4 . Thus during rotation 
fault normals move from point 2 to point 3. In time, with increasing stress, 
the fault normal will continue to rotate. When point 4 is reached (Fig. 8 c), 
the Coulomb fracture criterion will come into play — the fault will be locked at 
position < 23 , and a new set of faults will be created (point 5 in Fig. 8 c). The 
path followed from point 2 to point 4 is the rotation path for Nur’s 2D model. 

Unfortunately, a 2 D model is limited to analysis of the special case where 
cr 2 lies within the preexisting fault set. To model faults rotating more generally 
within the stress field, we must study the 3D case. 

3. 3D rotation paths 

Many rotation paths are possible in 3 D, but we will consider only three 
limiting cases, shown in Fig. 9. As described earlier, faults slip and rotate in 
the BR model, as long as the stress limits and the Coulomb sliding criterion 
are satisfied. By assuming a constant (f> stress history, we can plot the rotation 
path of a fault set on Mohr circles normalized by < 7 i. Thus, rotation paths 
plot as curved lines, reflecting the variation, with fault rotation, in the absolute 
magnitudes of the principal stresses. 

For a given stress regime, we would like to determine the style of faulting 
exhibited at every point along the rotation path of a fault set. This style can 
be characterized by the rake, or the direction of slip of the fault’s hanging wall 
along the footwall, as shown in Fig. 10. Because the fault plane itself is also 
rotating, we must track this change as defined by fault plane strike and dip, or 
more simply by the orientation of the fault normal in 3D space. 

In the following section, we present the results of 3 D BR modeling. We have 
assumed stationary stress orientations and constant <j> during rotation. Sepa- 
rately, we discuss modeled changes in faulting style (rake) and attitude (strike 
and dip) for the three specific fault sets, as functions of both the stress ratio <j) 
and the tectonic stress regime. 

3D BLOCK ROTATION: RESULTS 

Rotation paths for the three limiting cases (Fig. 9) are summarized in Fig. 
11 and Fig. 12. In each figure, results are shown for the three stress regimes, at 
low, intermediate and high <f> values. Fig. 11 displays changes in faulting style as 
a function of rotation in its nine rotation /rake plots. The columns represent the 
three <f> cases, and the rows represent the three tectonic regimes. The faulting 
style along the rotation path, that is the rake of a fault set, is determined by 
the fault set orientation and the <f> value: the higher the <f> value the more stable 
the rake becomes during rotation. This is illustrated by the Mohr circle plots 
(bottom row in Fig. 11 ): as the the o 2 — <73 circle grows, the 0 \ — a 2 circle 
shrinks, and rotation paths stop farther away from the a 1 — o 2 circle. 

Fig. 12 displays the fault set attitude for the same rotation paths shown in 
Fig. 11 . The columns still represent the three <j> cases, but Fig. 12a, b,c represent 
the normal, reverse and strike-slip regime seperately. 3D motion is shown by 
plotting rake vs. strike, rake vs. dip and strike vs. dip along the rotation path. 
Rotation paths are plotted on Wulff-projections at the bottom of Fig. 12 , to help 
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visualize the changing orientation of fault sets during rotation in a stationary, 
<{> = constant stress field. Notice, again, the stabilizing effect of high <f> values 
displayed by the shape of the rotation path. The dominant <J\ at <p = 0.1, drives 
all three paths very close to its direction. Some curvature of the paths can be 
noticed at <f) = 0.5, while straight paths characterize the <f) — 0.9 case — where 
<j\ and (72 are nearly equal in magnitude. 

Let us consider each rotation path in more detail. The initial orientation of 
pole 3 in Fig. 11 and 12, for example, closely contains the <t 2 direction in the 
planes of its fault set. This makes it equivalent to the 2D block rotation case — the 
intermediate stress <t 2 plays no part in this pole’s rotation path. The rake (Fig. 
11) remains unchanged dining rotation under any value of <f>: normal regimes 
produce normal slip. Likewise, reverse and strike-slip regimes yield reverse and 
strike-slip motion, respectively. Thus, our assumption that one principal stress 
axis is vertical has precluded oblique slip in the 2D case of pole 3. 

Pole 1 and pole 2, however, are influenced by <j 2 . As a result, their faulting 
style changes under rotation — sometimes dramatically. Surprisingly complex 
patterns of nonlinear rotation paths are outlined for these fault sets in Fig. 11 
and 12. Note that these changes take place under fixed <f> values and stationary 
principal stress directions. 

Changing faulting style — Rake 

As mentioned earlier, higher <f> values have a stabilizing effect. Thus, rakes 
for the rotation paths of pole 1 and 2 in Fig. 11 are nearly constant during 
rotation in all three stress regimes. At lower values of <j > , however, the model 
yields surprising rake histories for the rotation paths of these poles. In the 
normal stress regime, strike-slip fault sets may rotate into pure normal fault 
sets. In the reverse stress regime oblique-slip fault sets with a large component 
of strike-slip may rotate into a reverse slip fault set (pole 2) or a strike-slip fault 
set (pole 1). 

The strike-slip stress regime displays the most dynamic results of the 3D BR 
model.The results suggest that it is possible for pure normal fault sets (pole 
1) to rotate into a strike-slip one, and subsequently into a reverse-oblique slip 
fault set — all within a strike-slip regime! Even at intermediate <f> values, dramatic 
changes take place in the rakes of pole 1 and 2, although they converge to strike- 
slip motion upon rotation. High </> values limit rake changes during rotation, as 
always, but pole 2 does tend towards a more strike-slip motion. 

The 3D BR model predicts rotations of up to 75° (each rotation step corre- 
sponds to 5° of rotation) for the most poorly oriented fault sets. By considering 
these limiting cases, we have modeled rotations much greater than those of Nur 
et al. (1986), and we have also found cases where the faulting style may reverse 
duuring the rotation path of a single fault set. 

Changing fault orientation — Strike and Dip 

Fig. 12 has been designed to express 3D rotation about oblique axes as as a 
function of dip, strike and rake — to better visualize fault rotation in 3D. Again 
the stabilizing effect of higher d> values is observed in all three stress regimes 
(Fig. 12a, b,c). 

Results of the model for the normal stress regime are detailed in Fig. 12a. 
Considering the low <f> column on the left, one can see how all three paths change 
mostly in dip. This horizontal axis of rotation is what one might typically expect 
for a normal stress regime, and the 3D BR model concurs. Intermediate <f> values 
suggest greater changes in strike during rotation, while at high <f> values the 
model predicts — surprisingly so — that strike changes more than dip, particularly 
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for pole 1. Vertical-axis rotation, usually associated with strike-slip regimes, is 
therefore possible in normal stress regimes as well. 

The results for the reverse stress regime shown in Fig. 12b, predict rake 
changes throughout rotation at low <j> values. Poles 1 and 2 are expected to 
change more in strike than in dip. Again, we see how the 3D BR model pre- 
dicts vertical-axis rotation of fault sets — this time in a reverse stress regime. 
Increasing <f> values induce predominantly steepening of the fault planes during 
rotation, resulting in the more familiar horizontal-axis rotations. 

The results for the strike-slip stress regime shown in Fig. 12c, represent the 
most striking results of the 3D BR model. Interestingly, the low <j> case allows 
for decreasing dip during normal slip, followed by a steepening of the dip when 
the rake becomes reverse. In the strike-slip regime, the model predicts that most 
fault sets will follow an oblique axis of rotation. One notable exception, though, 
can be found in the <f> = 0.9 case, where pole 1 rotates unexpectedly about a 
horizontal axis in a strike-slip regime. 

To summarize the results of our model, three specific initial orientations of 
fault sets were studied in nine different combinations of (j> values (<f> = 0.1, 0.5, 0.9) 
and stress regimes (S v = <7j, ct 2 , cr 3 ). In each case the predictions of the model 
were presented using the rotation path. Along this path we analyzed behaviour 
(rake) and attitude (strike and dip) of the fault set as it rotated in a stationary, 
<t> = constant stress field. 

The most important conclusions of our results are: first, that for most cases 
the results predict an oblique axis of rotation, and second that during rotation, 
as faults change their orientation in the stress field, the style of faulting may 
change as well — sometimes dramatically. 

A direct consequence of these result is that paleomagnetically inferred ro- 
tations may not be directly related back to a specific tectonic stress regime. 
Indeed, rotations about vertical axes, while usually found in strike-slip stress 
regimes (Hornafius, 1986; Kamerling and Luyendyk, 1985; and others), have 
been documented in normal stress regimes as well (Li et al., 1990; Pavlides et 
al., 1988; Hudson and Geissman, 1987; Kissel et al. 1986; Brown and Golombek, 
1986; Jackson and McKenzie 1984; 1983). Paleomagnetism is an invaluable aid 
to decipher complex histories of rotation in regions of distributed deformation. 
Our 3D BR model provides a framework within which, paleomagnetic, struc- 
tural and stress data can be combined to better understand complex rotation 
histories. We analyze such a case in the following section, to demonstrate the 
application of the 3D BR model to actual complex tectonic problems. 

A 3D BLOCK ROTATION EXAMPLE: 

THE WESTERN TRANSVERSE RANGE 

The shear motion between the Pacific and the American plates in Southern 
California is distributed across a 200 km wide zone which consists of a complex 
array of block-faulted domains (Fig. 1; Luyendyk et al., 1985). In the first part 
of the paper, we showed (Fig. 7) how mixed styles of faulting observed today 
in this region can be explained in a stationary strike-slip stress regime. Now, 
we apply the results of the 3D BR model to the rotation history of one of these 
domains, namely the Western Transverse Ranges (WTR). 

1. Tectonic history 

The WTR domain is a region limited to the north by the Santa Ynez fault, 
to the south by the Malibu fault system, to the west by the Hosgri fault and 
to the east by the San Gabriel fault (see WTR in Fig. 1). The right-lateral 
shear between the Pacific and North American plates since Oligo-Miocene time 
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should have dominated the tectonic history of this domain. From the geological 
record alone, however, the tectonic history seems much more complex. This is 
documented by detailed stratigraphic and structural studies (Yerkes et al.,1981; 
Luyendyk et al., 1985; Yeats, 1987; Namson and Davis, 1988 and others). The 
following discussion is concerned with the most prominent deformational phases 
documented for this region — simplified for the purpose of this study. 

Phase I: in late Oligocene to early Miocene time, there was deposition of sed- 
iments in elongate normal-fault controlled basins which are presently oriented 
east-northeast (Terres and Lyuendyk, 1985). 

Phase II: a mixture of strike-slip and normal faulting characterized the exten- 
tional tectonics of Miocene times (Yeats, 1987). 

Phase III: a significant period of compression followed at the beginning of the 
Quaternary, overprinting all previous events. Structural evidence suggests (Yerkes 
and Lee, 1979a) that “at the present rates all the measured compressive defor- 
mation within the WTR could have occurred during the last 0.5 to 1 m.y.”. 
Namson and Davis (1988) describe this last phase of north-south convergence in 
detail. 

This sequence of phases is shown schematically in Fig. 13: a rake of 90°, 
representing the normal faulting period — phase I is followed by a rake of 0°, 
representing the strike slip faulting period — phase II, leading to a rake of —65°, 
representing the present day period of reverse oblique faulting — phase III. 

With the advent of paleomagnetic studies carried out in many domains across 
Southern California (Luyendyk et al., 1980; Terres and Luyendyk, 1985; Hor- 
nafius, 1986: Carter, 1987), it appears that block rotation has been the predom- 
inant mechanism of deformation in this region. Paleomagnetic interpretations 
indicate that since mid-Miocene time, the Transverse Ranges have experienced 
significant block rotations about a vertical axis — particularly in the WTR do- 
main, where up to 90° of rotation is estimated. 

According to Luyendyk and his co-workers, the WTR clockwise rotations 
were associated with left-lateral slip on a set of vertical, strike-slip faults origi- 
nally N-NE trending. The faults defined blocks about 100-200 kms long and 10 
kms wide that rotated as rigid bodies away from the direction of compression. 

The structural history for the WTR, as suggested by Terres and Luyendyk 
(1985), is explained by three separate set of faults formed under three different 
stress regimes (Fig. 13). A set of N-NE striking normal faults formed during 
the extensional phase (I), another set of strike slip faults, which rotated with 
the blocks, formed during the shearing phase (II), and yet another set of E-W 
striking high-angle reverse faults formed during the present-day compressional 
phase (III). 

While kinematic knowledge has improved our understanding of the WTR 
structural history, the story it tells is complex. By combining the frictional 
constraints of the 3D BR model with paleomagnetic, structural and geological 
data, we can now show how one set of faults, preexisting and rotating in a 
stationary strike slip stress field, can account for all three deformational phases. 

2. The 3D BR model for the WTR 

How can we use the results of Fig. 11 and 12 to help us understand the 
history of rotation of the WTR? Let us start with the simplest assumption: the 
stress directions remained stationary, in the present day orientation, throughout 
the post-Oligocene deformation history of the domain. From the first part of 
the paper, we know that it is possible for the WTR fault set to slip in a oblique 
reverse motion in a strike-slip regime with <Ti N20°E directed, if fit = 0.4 and 
<f) = 0.1. Consequently, we can expect the modeled rotation path that best 
describes the deformation of the WTR domain to be similar to that of pole 1 in 
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Fig. 11, for the case of the strike slip regime and <j> = 0.1, where the final value 
of the rake is oblique reverse. 

The results shown in Fig. 11 and 12 are for nj = 0.6. To model the faults of 
the WTR domain we must re-run the program with fif = 0.4. Fig. 14a shows 
the stress regime and the rotation path modeled for the WTR domain. Fig. 
14b shows how lowering the coefficient of frictional sliding allows faults to rotate 
further away from a x . The shape of the rotation path remains the same, but the 
limits within which a fault set may rotate change. In our case, faults may rotate 
until a steeper dip is achieved (compare strike and dip for the final orientation 
of pole 1 in Fig. 12c and Fig. 14b). 

But what does this rotation path tell us about previous fault behavior and 
previous geometries of faulting? Fig. 14b displays the structural history pre- 
dicted by the model for this rotation path. The rake starts at an initial value near 
pure normal slip, a time of “extension” similar to phase I. Upon fault and block 
rotation, the rake becomes pure strike-slip, a time of predominantly vertical axis 
rotation similar to phase II, before reaching the present day “compressive” pe- 
riod of reverse-oblique faulting — phase III. The corresponding values for the dip 
and the strike of the rotating fault set are also shown. As the strike of the fault 
set rotates away from <j\ (I — ► II in the plot), the dip flattens at first and then 
steepens as soon as the reverse slip field is reached (II — *■ III). 

3. Agreement with the observations 

Our goal was to demonstrate how deformation in the WTR could have oc- 
cured along a single set of faults, reactivated in middle Miocene time as block 
rotation began. 

In agreement with the observations (Fig. 13), Fig. 15 displays the salient 
features predicted by the 3D BR model for the structural history of the WTR 
fault set. The original orientation of the fault set in the WTR domain must 
have been NNE striking, dipping 55° — 60° to the West and slipping in a normal 
sense. Subsequently, the same fault set rotated away from the direction of 
maximum compression a x and became strike-slip in style with a dip between 40° 
and 50°. Finally, the fault set rotated into the E-W, high angle, reverse-oblique 
faults observed today. The results of our modeling show that along the rotation 
path, a single fault set went through three different faulting styles during its 75° 
clockwise rotation, while stress directions remained stationary with S v = <72, a \ 
oriented N20°E, and <f> = constant at a value of 0.1. 

These simplifying assumptions are not necessary. If we most accurately mod- 
eled the tectonic history of the WTR domain, we should have accounted for 
changes in stress directions that presumably took place between the Miocene 
and the Present. Evidence of this is found in studies of plate motion (Cox and 
Engebretson, 1985) that indicate a clockwise change in relative plate motion 
~ 5 My ago. Incorporating this change would introduce an additional degree of 
freedom and allow a less restrained solution. The original orientation of faults 
could be even more N-S striking, and the amount of total rotation that could be 
achieved by one single set of faults would increase by an amount comparable to 
the change in plate motion. 

The second constraining assumption, that of a constant <f> stress history, 
mantains a strike-slip stress regime during the rotation history of the faults. 
This need not be the case. A change of stress regime might occur dining the 
history of rotation of a fault set, and it could reverse the sense of rotation of 
blocks and faults or allow them to rotate in another direction. 

Nevertheless, it is encouraging that with our simple 3D BR model and sim- 
ple assumptions, we have been able to predict a sequence of structural phases 
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consistent with the known structural history of the WTR. This example demon- 
strates that one need not always invoke complex regional and local changes in 
the stress regime or erratic changes in plate motions to account for alternate 
periods of compression and extension. The geometry of each domain of fault 
sets determines the style of faulting that will occur there. The regional stress 
regime may remain stationary. 

CONCLUSIONS 

We have presented a generalized model of distributed deformation and block 
rotation that emphasizes faulting in 3 dimensions. The need for this model is 
dictated by the following observations: 

fl) Reactivated faults are usually poorly oriented relative to the stress field. 

(2) Obliquely slipping faults are found in many domains that have undergone 
rotation. 

(3) These same faults change their faulting style through time. 

We propose that mechanical constraints of friction and strength combined with 
the 3D Mohr circle and the kinematics of block rotation can explain these ob- 
servations. 

1. Block rotation: the cause of poorly oriented faults 

According to the existing 2D BR model, rotation of blocks and faults are 
caused by contemporaneous slip on fault sets in adjacent domains. Based on 
fracture mechanics and friction criteria, the orientation of faults sets with respect 
to the direction of maximum compression o\ determines the sense of fault slip 
and thus the sense of rotation of the blocks in each domain. Consequently, faults 
that are initially in optimal orientations must rotate away from this optimal 
orientation in domains undergoing rotation. In this paper, we generalized these 
concepts into a 3D BR model. 

With the 3D BR model, we are able to determine the range of orientations 
for faults that may slip as well as their style of faulting (i.e. rake). We can also 
determine the maximum amount of rotation blocks and faults may experience 
before a new set of faults may take over the deformation. 

The results of our modeling show that the value of the 4> parameter plays a 
key role in identifying those faults that may be reactivated and in determining 
their style of faulting. At low <j> values, where <r 2 is close in magnitude to 
cr 3 , faults may be reactivated even if very poorly oriented. Thus, in different 
domains at the same point in time or within the same domain through time, 
contrasting styles of faulting can occur. In a normal stress regime, both pure 
strike-slip and pure normal faulting styles can occur. In a reverse stress regime, 
both pure strike-slip and pure reverse faulting styles can occur. In a strike-slip 
stress regime, pure strike-slip, pure normal and pure reverse faulting styles can 
all occur. At high <j> values, where cr 2 is closer to cr x in magnitude, faulting styles 
are limited to strike-slip and normal in both the normal and strike-slip stress 
regimes and only to reverse in the reverse stress regime. 

When applied to the distributed deformation of Southern California, our 
model can explain the diverse faulting styles observed there by assuming a simple 
regional strike-slip stress regime with 0 \ oriented N20°E, a low (f> = 0.1, and a 
coefficient of sliding friction of 0.4. 

2. Change of faulting style with rotation 

The importance of the <f> parameter is also clearly demonstrated by the change 
in rake predicted by the 3D BR model as faults rotate in a stationary stress 
field. The results show that rake varies with rotation at low <f> values, but 
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remains relatively constant at high values. By analyzing dip and strike of faults 
throughout rotation, we found surprising results. At low <j) values in the reverse 
regime, originally poorly oriented faults may rotate about a vertical axis. On 
the other hand, at high <f> values, poorly oriented faults may rotate unexpectedly 
about a vertical axis in the normal stress regime, and about a horizontal axis in 
the strike-slip stress regime. Therefore, for the general case of originally poorly 
oriented sets of faults, oblique axes of rotations are predicted by our model for 
all stress regimes. 

We applied the 3D BR model to the complex history of faulting as doc- 
umented in the Western Transverse Ranges of Southern California. The best 
model results indicate that for the WTR domain, a single set of originally NNE 
striking faults could have accommodated the deformation. Assuming a low <f> 
strike-slip stress regime, a stationary N20°E directed maximum compressive 
stress, and a value of 0.4 for the coefficient of sliding friction, these faults were 
reactivated as normal faults at the beginning of the rotation period. Subse- 
quently, the same faults became more NE striking, left-lateral strike-slip faults 
until they finally rotated into the present day “compressive” period of E-W 
striking reverse faults. 

In spite of its simplifying assumptions, the our model provides clues to ear- 
lier geometries of faults for the WTR that are consistent with the geological, 
geophysical and paleomagnetic records. 

We believe that block rotation is a fundamental process in regions of dis- 
tributed deformation where poorly oriented faults are widespread. When ap- 
plied to actual geologic situations, the model provides a powerful, yet simple 
tool to interpret the complex faulting histories and complex fault geometries 
that characterize regions of distributed deformation. 
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APPENDIX 


1 . Stress limits 

A fault can slip and rotate when the Coulomb criterion for sliding is satisfied: 

T f = C/ + flf(T n ( 1 ) 

where r is the shear stress on the fault plane, a n is the normal stress across it, 
Cj is the cohesion and Hj is the coefficient of friction of the fault. 

Prior to slip and rotation, it is necessary to check that the differential stress 
does not exceed the Coulomb fracture line: 


T 0 — C 0 + t L O a n (2) 

where C 0 is the cohesion and fx 0 is the coefficient of friction of the intact crust. 
This line defines the upper limit for the magnitude of the differential stress (R) 
in the earth’s intact crust : 

R — X < Rmax ~ J= -[Co + 1 + CT 3 )] ( 3 ) 

2 Vl + Z^o 2 

Furthermore, it is necessary to ensure that: 


<r 3 > 0 (4) 

because the crust cannot support tension. 

2. Stress model and stress paths 

The normal (<r n ) and shear (r) stresses across a plane are computed according 
to the following relations (Jaeger and Cook, 1969): 

<7 n OjCj + 0i 2 0 2 ”1“ ^3*^3 (5) 

r 2 = a\cr\ + a 2 2 a 2 2 + alal - <r 2 (6) 

where ari , a 2 , 03 are the direction cosines of the pole (fault normal) with respect 
to the principal stress axes, <r l5 <r 2 , <r 3 . 

As principal stress magnitudes change to satisfy equation (1), the initial value 
(^n, t) on a given fault plane will change to a new value (a' n , r'). The sequence of 
(<7' n , t' ) described by the pole is computed with «i, a 2 , a 3 = constant , because 
during these stress changes and before any rotation event occurs, the fault plane 
itself remains fixed in stress space. 

This sequence of (<T^,r / ) will be referred to as the stress path of a fault. It 
depends on the stress history chosen and on the initial conditions (cr„, r, 

03, ot 1, a 2 , 03)’ 

As explained in the text, this paper considers three stress regimes in which 
one of the principal stress axis is always assumed vertical (S„) and the stresses 
are assumed to change according to a <j> = (o 2 - cr' 3 )({a[ - (t' z ) = constant model. 
For each stress regime, the stress path of a fault, is computed as follows: 

a. NORMAL STRESS REGIME (S v = a t ) 

o\ — a\— constant ( 7 ) 
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where 


o' 2 = °z ~ r^Aa 
<73 = °3- Act, 


( 8 ) 

( 9 ) 


r N 


a l ~ a 2 


<J\ — ct 3 


b. REVERSE STRESS REGIME ( S v = ct 3 ) 

cr[ = cr 1 + A (7 


where 


^2 = ^2 + FrAcT 
a 3 = a 3 = constant 


02 ~ 03 
r R = 

V 1 — 03 


( 10 ) 

( 11 ) 

( 12 ) 

( 13 ) 

(14) 


c. STRIKE-SLIP STRESS REGIME ( S v = cr 2 ) 

a\ = era -f Act (15) 


^2 ~ <J 2 — consianf 


(16) 


03 = 03 ~ r S A <7 

where 

(72 — 0*3 

r 5 = 

<Ji — cr 2 


( 17 ) 

(18) 


Substituting the values for o\ into equations (5) and (6), the stress paths 
that result from the above models are given by the following equations with Act 
as a parameter: 


cr' n — o n + aAo 

(19) 

= t 2 + bAo + c(Act) 2 

(20) 
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where <r 1? cr 2 , <j 3 , r, cr n are known initial values and Qj , a 2 , 03 are the constant 
direction cosines of the fault normal. The coefficients a, b and c assume different 
values depending on the stress regime considered (see Table Al). 



Coefficients 

Stress Regime 

a 

b 

c 

normal 

~(<*3 + r N<* 2) 

-2[r N a\(a 2 - a n ) + a\(a z - cr n )] 

(1 - r N ) 2 ala 2 + (a 2 + r 2 N a 2 2 )a 2 

reverse 

a 2 + r R a\ 

2 ( Q i( cr i - ^n) + r R al(a 2 - <r n )] 

(1 - r R ) 2 a 2 a 2 -f (a? + r 2 R al)a\ 

strike-slip 

a i ~ r S a 3 

2[ a i( cr i - a n) ~ r s al(a 3 - <r n )] 

a? + r 2 a 2 - a 2 


Table Al. a, b and c are coefficients required to calculate the stress path of a fault. 
They depend on the orientation of the faults ( ot\ , a 2 , 013), the assumed stress 
model (Vjv, r R , rs) and the initial values of <j\, a 2 and <73. 


Replacing A a from (19) into (20), we get an explicit equation for the stress 
path 


a 2 r n = Aal + Ba' n + C, (21) 

where 

A = c (22) 

B = ab- 2ccr n (23) 

C = a 2 r 2 - abcr n + ca 2 n (24) 

By combining equation(21) and equation(l), we obtain the coordinates cr' n ,r' 
of the intersection between the pole’s stress path and the Coulomb sliding fine. 
There the fault plane can slip and rotate. 

3. Rotation of fault planes 

Consider a normal to a fault plane with unitary vector direction cosines 
a ii a 2 5 <*3- Once the shear stress r across the fault plane satisfies equation (1), 
the fault plane will slip and rotate in the direction of the resolved shear stress, 

represented by the unitary vector f. 1 The nature of this rotation is better 

understood in Fig. Al, where the horizontal plane represents the fault before 
rotation and r is its unit normal vector (same direction as a n ). 

Given the coordinates o^o^, 03 of the normal to the fault r and the coor- 
dinates of the total stress vector (a) across it 

[£]i = <*1^1 (25) 

[a] 2 = a 2 <7 2 (26) 


1 A tilde will denote a unitary vector 
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[a] 3 = a 3 cr 3 , (27) 

we can compute the coordinates of vectors h, r and f . In order to compute h 
we first compute vector n = fx a with coordinates: 

nj = a 2 a 3 (a 3 - a 2 ) (28) 

n 2 = aia 3 (ai - <r 3 ) (29) 

n z = aia 2 (<7 2 - <7 X ) (30) 

Vector n is then simply n/\n\. Similarly, ?=n x r and has coordinates: 

T X = ai {al(<T! - <r 3 ) - a|(a 2 - crj)} (31) 

r 2 = a 2 {aj(a 2 - a x ) - a|(<r 3 - <r 2 )} (32) 

r 3 = a 3 {Q2(^3 - 02 ) - a?(<7-i - 0 - 3 )} (33) 

The direction cosines of the rotated plane will be the coordinates of the vector 

r'. Therefore, they must satisfy the following 3 equations: 

r' • n = 0 => Oi[n\ + a' 2 n 2 + Q: 3 n 3 = 0 (34) 

r' ■ r = cos 6p aiai + 0 : 20:2 + <* 3 <*3 = cos<5p (35) 

r'-r = r'-r/lrl = sin^p =$> ajri + a 2 r 2 + <* 3 73 = sin^p • |r| (36) 


After rotation has occured, principal stresses will drop and Mohr circles should 
“shrink”. The newly computed values of the rotated fault plane direction cosines 
(aj,a 2 ,o: 3 ) define a new stress path, along which the pole must move as the 
stresses change (after each rotation we can choose a different model for stress 
changes). No direct relationship exists yet between the amount of stress drop 
and the amount of rotation. Nonetheless, whatever the stress drop, the new pole 
must move along its stress path. Since the rotation path is fully defined by the 
intersection between the stress path and the sliding line, it is not necessary to 
know the stress drop (that is, how far the pole will move along its new path) to 
estimate the next rotation event. 
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Fig. 1 Structural domains in Southern California. Sub-parallel 
sets of faults (Quaternary in Age) of different orientations define structural 
domains: Mojave (MOJ) and East Mojave (EMOJ), West (WTR), Cen- 
tral ( CTR ) and East Transverse Ranges ( ETR '.). Domains are separated by 
throughgoing boundary faults such as the San Andreas, Santa Ynez, and Gar- 
lock faults (base map adapted from Jennings, 1975). Representative fault 
plane solutions (W, E and M) for three domains and for the San Andreas 
fault (A) are also shown. 


2b 





Fig. 2 Kinematics of block rotation. Blocks (books) of different orien- 
tations abut against boundary faults (shelf and book-holders) . The orientation 
of the blocks with respect to the direction of maximum compression (in this 
case gravity) controls the sense of slip and therefore the sense of rotation of 
blocks. 



Fig. 3 Optimally oriented fault. According to traditional friction 
models (Zoback and Healy, 1984), active faults are generally optimally orientd 
in the principal stress axis. Therefore, when plotted on the Mohr circle , 
optimally oriented faults should be represented by a <7\ — <73 circle tangent to 
the sliding line , with the intermediate stress 0*2 contained in the plane of the 
fault. 




Fig. 4 Mechanics of block rotation. The 2D Mohr circles show values 
of <7 n and r required for the fault to slip , when its fault normal is at point 1 , 
and when the fault rotated to a point 2 , closer to the direction. The lines 
represent the Coulomb criterion for sliding along preexisting faults and the 
Coulomb criterion for fracturing an intact rock mass. As illustrated in the 
sketch to the right, faults rotate away from cf\. Therefore, normals to faults 
(poles to fault planes) will rotate towards cr 1; on a Mohr circle. 
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Fig. 5a The SD Mohr circle is a S-dimensional construction on a 2- 
dimensional a n — r plane. Any point in the shaded region represents the 
oHentation of a fault plane in space. For example, aj and <*3 are the direction 
cosines of a pole (P) relative to the <j\ and <73 directions respectively. Three 
specific cases are shown: pole (Pi) plotting on the (J\ — <73 circle with the 
<t 2 direction contained in its fault plane, identical to the 2D case; pole (P x ) 
plotting on the a 2 — <r 3 circle, with the 0\ direction in its fault plane and pole 
(P3) plotting on the U\ — a 2 circle, with the 03 direction in the fault plane. 



Fig. 5b Poles P, P X) P 2 and P 3 of Fig. 5a are plotted in a lower hemi- 
sphere Wulff stereo projection. This is another way of plotting the orientation 
of poles to faults in a SD space. 
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Fig. 6 Mohr circle, friction criteria and faulting styles. Different 
faulting styles (N = pure normal, R = pure reverse, S = pure strike- 
slip) may occur in a region of distributed deformation. It depends on the 
orientation of the fault sets, the tectonic stress regime and the stress ratio 
(f> as well as on the material parameters. Here we consider three sets of 
faults subjected to the three tectonic stress regimes ( assuming 1.0 and .6 for 
the fracture and sliding coefficeints respectively): set 2 which contains the o 2 
direction in its plane 1 the 2D case), set 1 which contains the <j\ direction in 
its plane and, set 3 which contains the a 3 direction in its plane. 
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Fig. 7 Representative fault planes shown in Fig . 1 are plotted on a 3D 
Mohr circle (normalized by o\). A slope of 0.4 for the sliding line , a value 
of <f> = 0.1 for the stress ratio and a regional strike slip stress regime with 
G\ N20°E directed are required to slip faults in the different domains. No- 
tice that the San Andreas fault (SAF), a major throughgoing boundary fault, 
cannot slip in accordance with friction criteria assumed in this model . 
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Fig. 8a. Assume a set of faults exists in a domain with an orientation aj 
to the <j\ direction. The direction is contained in the fault plane. Stress 
magnitudes must change (1 — » 2) to allow sliding on the fault set. Assume 
that this change takes place by keeping the mean stress cr 0 constant. 



Fig. 8b. The Rotation Path; At point 2 there is sufficient shear stress 
to overcome friction. The faults of a set slip and rotate. On the Mohr 
circle poles to faults are plotted (see Fig. 4). Thus, as rotation proceeds, 
stress magnitudes change and the poles plot closer and closer to o x tracing a 
rotation path from 2 to 3 along the sliding line. 
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Fig. 8c At point 4 rotation along this set of faults ceases once stress 
magnitudes have reached the fracture line . Therefore a new , more optimally 
oriented set of faults forms, point 5 } and the old set remains locked into 
an orientation a 3 relative to <j\ . In 2D, only fault orientations with the a 2 
direction in their plane can be modeled. After each rotation step , faults plot 
on the < 7 i — circle. This results in only one rotation path and a constant 
fault behaviour. 



Fig. 9. In 3D there are many rotation paths . In this paper we consider 
three specific paths. They all start close to the <?$ direction but rotation path 
for pole 1 starts close to the a 2 — 03 circle, rotation path for pole 3 starts 
close to the — cr 3 circle, rotation path for pole 2 starts at an intermediate 
position between the two. Rotation paths shown in this figure and Fig. 11 
are plotted on graphs normalized by the magnitude of <j\. The stress ratio <fi 
remains constant throughout the rotation history of a fault set. Thus rotation 
paths will appear as curved lines. 
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Fig. 10: The orientation of a fault set is given by the strike , or the 
azimuth of the fault with respect to a reference direction , and the dip, or 
the angle between the fault plane and the horizontal plane (0° = horizontal; 
90 ° = vertical). The style of faulting is determined by the rake, or the angle 
between the direction of fault slip and the horizontal (+ 90 ° = pure normal; 
0 ° = pure strike — slip; — 90 ° = pure reverse). 
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Fig. 11 Cross-plot of faulting styles for different stress regimes and <f> 
ratios. The faulting style evolves along various rotation paths in plots of rake 
vs. rotation. Strike-slip fault styles (S ) are at 0°, normal (N ) styles are at 
90°, and reverse (R) styles are at —90°. Fig. 9’s case 1 plots with a solid line, 
case 2 with dashes, and case 3 with dots . Variations in rake with rotation 
are sensitive to cr 2 . Low <f> values change the modeled rake dramatically while 
high values leave it largely unchanged along the rotation paths. 
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Normal 



Fig. 12a Rake versus Dip, Rake versus Strike and Dip versus Strike are 
plotted for the three rotation paths (1, 2, 3) at different values of the stress 
ratio <j>. The arrows indicate the sense of rotation. Wulff stereo-projections 
for the three rotation paths are shown, for reference, at the bottom of each 
figure. 

Normal stress regime, S v = <j\. The flattening of the dip decreases with 
increasing stress ratio. The rake changes from strike-slip to normal at low 
stress ratios (f> = 0.1) but remains constant at the initial value at high stress 
ratios. This is accompanied by an increase in vertical axis rotation with in- 
creasing stress ratio (see rotation paths 1 and 2): the strike is relatively 
constant at low stress ratios but changes considerably at high ratios (a$ is at 
strike= 0°). 
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Reverse 



Fig. 12b Reverse stress regime, S v = <r 3 . The dip steepens through- 
out rotation for all fault planes at all stress ratios. The rake changes from 
oblique-strike-slip to oblique-reverse-slip (rotation path 1) or to strike-slip 
(rotation path 2) at low stress ratios but at intermediate and high stress ra- 
tios the rake is reverse along all rotation paths. This is accompanied by a 
vertical axis rotation at low 4> values (see rotation paths 1 and 2) where the 
strike changes considerably. At higher stress ratios, rotations occur predom- 
inantly about a horizontal axis (&2 is at strike= 0°). 
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Fig. 12c Strike slip stress regime, S v = cr 2 . As the rake goes from 
normal , through strike-slip and then to reverse, the dip of the fault plane 
flattens at first and then steepens at low stress ratios (<f> = O.lJ. At higher 
stress ratios rotation stops before the reverse field is reached. The strike 
change is greater than the dip change, producing predominant vertical axis 
rotations in most cases (o\ is at strike= 0). An exception to this is rotation 
path for pole 1 at <f> = 0.9: the strike does not change, only the dip and 
therefore a purely horiiontsl <ms rolahon can occur in a strike-slip stress 
regime. 
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Rake 





Phase 


Phase 


Phase 


Fig. 13 Observed and inferred rake, dip and strike of the WTR fault set 
since Miocene times. Based on structural, geological and paleomagnetic data 
as discussed in the text. o\ si assumed at N20°E and is at strike= 0° on the 

plot. 



Fig. If Modeling 3D fault rotation in the West Transverse Ranges. 
(a) The rotation path predicted by the SD BR model, assuming, as in Fig. 

7, a strike slip stress regime with the maximum principal stress direction a x 
N20°E, a coefficient of friction pj = 0.4, and a stress history with a constant 
stress ratio <f> = 0.1. (b) According to the 3D BR model, a single set of ro- 
tating faults may have accommodated the deformation of the WTR domain. 
The model predicts, in accordance with observations summarized in Fig. 13, 
three different phases of faulting styles since rotation began: normal-oblique 
faulting (I) — > strike-slip faulting (II) — ► high angle reverse-oblique faulting • 
(III). 
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Fig. 145 The evolution of the rake, the dip, and the strike of the faults 
of the WTR through time. Each phase may be better visualized with idealized 
block diagrams and focal plane solutions (lower hemisphere projections). The 
sketches in the right column show the dominant axis of rotation, where S v 
refers to the principal vertical stress axis and the strike of the fault is shown. 
When the strike of the fault is almost parallel to C\ — phase I — the faults move 
in a normal-oblique (Vake= -\-70) sense. When none of the principal stresses 
are parallel to the strike of the faults — phase II — the faults move in a strike- 
slip sense frake= 0). Finally, when the strike of the fault becomes parallel 
with the cr 3 direction — phase III — the faults move in a reverse-oblique sense 
frake= -65). 


37 







Fig. Al: Rotation of Fault Planes. The horizontal plane, with its 
normal vector r, represents the initial orientation of the fault. During rota- 
tion, the normal r will rotate through an angle 8p within the plane (normal 
h) that contains <j n , r, a , t, moving to a new position V , defined by direction 
cosin es a j , a' 2 , 03 . 
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